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Abstract 

We complete our study of non-Abelian gauge theories in the framework of Epstein-Glaser approach 
to renormalization theory including in the model an arbitrary number of Dirac Fermions. We 
consider the consistency of the model up to the third order of the perturbation theory. In the 
second order we obtain pure group theoretical relations expressing a representation property of 
the numerical coefficients appearing in the left and right handed components of the interaction 
Lagrangian. In the third order of the perturbation theory we obtain the the condition of cancellation 
of the axial anomaly. 
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1 Introduction 



In some preceding papers [16], [17] we have extended results of Aste, Diitsch and Scharf ||, [13], 
concerning the uniqueness of the non-Abelian gauge theory describing the consistent interaction of 
Bosons of spin 1. It appeared that the gauge invariance principle is a natural consequence of the 
description of spin-one particles in a factor Hilbert space: gauge invariance expresses the possibility 
of factorising the S-matrix to the physical space, which is usually constructed using the existence 
of a supercharge Q according to the cohomological-type formula: H p h ys = Ker(Q)/Im(Q). The 
obstructions to such a factorization process are the well-known anomalies. The case when the spin- 
one Bosons of non-null mass are admitted in the game was studied in |L|] , for the concrete case of 
the electro-weak interaction i.e. when the gauge group is exactly SU(2) x U(l). 



In [17] we have analysed the same problem considering that the spin-one Bosons can have non-null 



masses and we did not impose any restriction on their number and masses and we did not took into 



account the matter fields. Similar results have been obtained in \24\- We have obtained only from the 
condition of absence of the anomaly up to the second order the existence of a Lie algebra g and the 
existence of a representation of this Lie algebra pertaining to the Higgs fields. 

In this paper, we consider the effect of including Dirac Fermions. In this way we are able to 
investigate a truly realistic model of gauge interactions of elementary particle and, in particular, to 
see what are the restrictions on such a model determined by the cancellation of all anomalies. The 
main results are the following ones. 

(A) The cancellation of the anomaly in the second order of the perturbation theory brings new 
relations on the numerical coefficients of the left and right handed components of the interaction 
Lagrangian. More precisely, new group theoretical property appear: 

(i) The coefficients of the vectorial and pseudo-vectorial couplings can be organised as two repre- 
sentations of the gauge algebra: t+ and t~ with a,b, . . . = 1, . . . , r group indices; the usual notations 
are and t^. 

(ii) The coefficients of the scalar and pseudo-scalar couplings can be organised as some tensor 
operators. 



Some of these relations have been obtained from different considerations in [26], ||. 
(B) The cancellation of the anomaly in the third order of the perturbation theory gives, essentially, 
the usual condition of cancellation of the axial anomaly: 

A abc = Tr (t+{t+ t+}) - Tr (t~{^,C}) • (1-0.1) 
This is the expression of the Adler-Bardeen-Bell-Jackiw anomaly Q, §, @, @, @, f20|, 



Some conditions on the couplings of the Higgs fields appear if one imposes the additional require- 
ment that no finite renormalizations of degree greater than 4 are allowed. This conditions gives the 



usual expression for the Higgs potential 13], ||] for the case of the Standard Model (SM). In [24] and 
[14] the analysis of the pure Boson sector is performed somewhat differently. We will comment on 
that in the next Section. 

The structure of the paper is the following one. In the next Section we define the model and 
construct the interaction Lagrangian including Dirac Fermions also. Then in Section || we outline the 
general setting for the study of the renormalization theory, the general structure of Ward identities 
and some facts about distribution splitting. In Section || we construct the S-matrix up to the second 
order of the perturbation theory. For the case without matter fields we reobtain the results of [17|. 



Then we consider the coupling of Yang-Mills fields with Dirac Fermions and, as anticipated above, 
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we obtain some interesting group-theoretical relations; their complete analysis could fix very severely 
the possible generalisations of the standard model. We also analyse the conservation of the BRST 
current in the second order of the perturbation theory. In Section |5] we go to the third order of the 
perturbation theory. We investigate the Dirac Fermionic sector and we get the new conditions on the 
Fermionic representations (a) or (b) from above. In the end we particularise the formalism for the 
case of the standard model with one generation of Dirac particles. In the last Section we comment on 
the possibility of removing the anomalies in all orders of perturbation theory. 

For the sake of clarity of the rather long and intricate analysis we adopt the mathematical definition 
- theorem style of presenting various assertions and computations. 
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2 General Description of the Vector Bosons 



2.1 Massive Yang-Mills Fields 



In |[(| and JT^] we have started from the following two facts: 

(1) a system of free zero-mass vector Bosons can be described in a Hilbert space generated from 
the vacuum f2 by applying the free fields A^, u, u of zero mass and a factorization procedure induced 
by a supercharge operator; 

(2) a system of free vector Bosons of mass m > can be described in a Hilbert space generated 
from the vacuum $7 by applying the free fields A^, u, u, of mass m and a factorization procedure 
induced by a supercharge operator. 

Here A^ is a Boson vector field, u and u are scalar Fermi fields and $ is a scalar Boson fields; the 
fields u, u, $ are usually called ghost fields. 

For the Yang-Mills model we combine somehow these two cases. We consider the auxiliary Hilbert 
space HyM generated from the vacuum Q by applying the free fields A afl , u a , u a , & a a = 1, ... , r 
where the first one has vector transformation properties with respect to the Poincare group and the 
others are scalars. In other words, every vector field has three scalar partners. Also u a , u a a = 
1, . . . , r are Fermion and A^, & a a = 1, . . . , r are Boson fields. 

We have two distinct possibilities for distinct indices a: 

(I) Fields of type I correspond to an index a such that the vector field A% has non-zero mass m a . 
In this case we suppose that all the other scalar partners fields u a , u a , <& a have the same mass m a . 

(II) Fields of type II correspond to an index a such that the vector field A% has zero mass. In this 
case we suppose that the scalar partners fields u a , u a also have the zero mass but the scalar field <3? a 
can have a non-zero mass: > 0. It is convenient to use the compact notation 



m a for m a ^ 
for m„ = 



(2.1.1) 



Then the following following equations of motion describe the preceding construction: 

(□ + ml)u a (x) = 0, (□ + m 2 a )u a {x) = 0, (□ + (m* a ) 2 )<S> a (x) = 0, a=l,...,r. (2.1.2) 

We also postulate the following canonical (anti)commutation relations: 

[A afl (x),A b „(y)] = -5 ab g^ v D ma (x - y) x 1, 

{u a (x),u b {y)} = S ab D ma (x-y) x 1, [$ a (x),$ h (y)] = 5 ab D m * (x - y) x 1; (2.1.3) 

all other (anti)commutators are null. 

In this Hilbert space we suppose given a sesquilinear form < -, • > such that: 

A afJ ,(x)^ = A afJ/ (x), U a (x)^ = U a (x), tt a (x) f = -U a (x), <5> a (x)^ = $a(x) . (2.1.4) 

The ghost degree is 1 for the fields u a , u a , a = 1, . . . , r and for the other fields. 
One can define the BRST supercharge Q by: 

{Q,u a } = {Q,u a } = -iip^ + ma^a) 

{Q,A^} = id"u a {Q,^ a } = im a u a , Va = l,...,r (2.1.5) 

and 

Qtt = 0. (2.1.6) 
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Then one can justify that the physical Hilbert space of the Yang-Mills system is a factor space 



H r YM = H = Ker{Q)/Ran(Q). (2.1.7) 

The sesquilinear form < , • > induces a bona fide scalar product on the Hilbert factor space. 
The factorization process leads to the following physical particle content of this model: 

• For m a > the fields Ag, u a , u a , <& a describe a particle of mass m a > and spin 1; this are 
the so-called heavy Bosons [jTfl| . 

• For m a = the fields Aa, u a , u a describe a particle of mass and helicity 1; the typical example 
is the photon 

• For m a = the fields $ a describe a scalar fields of mass m^', this are the so-called Higgs fields. 

This framework is sufficient for the study of the Standard Model (SM) of the electro-weak interac- 
tions: indeed one takes r = 4 and considers that there are three fields of type I and one field of type 
II. The scalar field appearing in the last case can be considered as the Higgs field. 

To include also quantum chromodynamics one must consider that there is a third case: 
(III) Fields of type III correspond to an index a such that the vector field Aa has zero mass, the 
scalar partners u a , u a also have zero mass but the scalar field 3> a is absent. 



In [24] and [14| the model is constructed somewhat differently: one eliminates the fields of type II 
and includes a number of supplementary scalar Bosonic fields tpi of masses mj > 0. In this framework 
one can consider for instance the very interesting Higgs-Kibble model in which there are no zero-mass 
particle, so the adiabatic limit probably exists. 

One can preserve the general framework with only two types of indices if we consider that in case 
II there are in fact three subcases (i.e three types of indices a for which m a = 0): 
(Ha) In this case A a ^, u a , u a , <3> a ^ 0; 
(lib) In this case <3? a = 0; 
(lie) In this case A afl , u a , u a = 0. 



One must modify appropriately the canonical (anti)commutation relations ( [2.1.3 ) to avoid con 



tradiction for some values of the indices. One has some freedom of notation: for instance, one can 
eliminate case (Ha) if one includes the first three fields fields in case (lib) and the last one in case 
(lie). The relations (gX| ) are not affected in this way. 

Let us consider the set of Wick monomials W constructed from the free fields Aa, u a , u a and & a 
for all indices a = 1, . . . , r; we define the BRST operator : W — > W as the (graded) commutator 
with the supercharge operator Q. Then one can prove easily that: 

a% = 0. (2.1.8) 

The class of observables on the factor space is defined as follows: an operator O : — > 
induces a well defined operator [O] on the factor space Ker{Q) / Im(Q) ~ T m if and only if it verifies: 
dQO\Ker(Q) = Because of the relation ( |2.1.8| ) not all operators verifying the condition ( |2.1| ) are 
interesting. In fact, the operators of the type dqO are inducing a null operator on the factor space; 
explicitly, we have: 

[d Q 0] = 0. (2.1.9) 

The canonical dimension w(H^) of certain Wick monomial is defined according to the usual pre- 
scription. By definition, a Wick polynomial is a sum of Wick monomials. 
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We will construct a perturbation theory d la Epstein-Glaser using this set of free fields and impos- 
ing the usual axioms of causality, unitarity and relativistic invariance on the chronological products 
T(xi, . . . , x n ). Moreover, we want that the result factorizes to the physical Hilbert space in the 
adiabatic limit. This amounts to 

= 0, Vn>l. (2.1.10) 

Ker{Q) 



lim e \ d Q / dxi---dx n g e (xi)'-g e (x n )T(x 1) ...,x n ) 

' J4^J xn 



If this condition if fulfilled, then the chronological and the antichronological products do factorise 
to the physical Hilbert space and they give a perturbation theory verifying causality, unitarity and 
relativistic invariance. 

One may raise at this point the rather serious objection that the adiabatic limit probably does 
not exists. One way to "cure" this problem is to replace the condition of factorisation ( 2.1 - 1C| ) by the 
"infinitesimal" version postulated in || - namely: 

n d 

d Q T(xx,...,x n ) =i^2-j— jI T l ti (x 1 ,...,x n ) (2.1.11) 

i=i X{ 

for some auxiliary chronological products T^(x\, . . . ,x n ), I = 1, ... , n which must be determined 
recurringly, together with the standard chronological products, and to construct the 5-matrix S(g) 
for a test function g, that is without performing the adiabatic limit g \ 1. 

However, this point of view is not without problems. Indeed, if one imposes (12.1.11|) instead of 



( 2.1.10 ), then the S-matrix so constructed will not factorize to the physical space Ker(Q)/Im(Q) 
which raises the question about its physical relevance. To this one must add the rather unpleasant 
fact that one abandons the consistency condition ( 2.1.10D which has a direct physical relevance (the 



possibility of constructing an 5-matrix in the physical space Ker(Q)/Im(Q)) for an independent 
postulate (]2.1,11 ). On the other hand, the rather close connection between (|2.1.10|) and ( [2.1,11 ) 



suggests that there must exists a common "cure" for both types of problems. That is, if one can find 
a reasonable solution of the adiabatic limit problem, then it is reasonable to conjecture that one will 
be able to strengthen the mathematical status of ( 2.1 ■ 10| ) and, eventually, prove its equivalence with 



( 2.1.11 ). In this case the consistency condition can be also written in the following form: 
dxi ■ ■■dx n g e (x 1 ) ■ ■ ■g € (x n )T(xi, ...,x n ) 



IQ 



= 0(e), Vn>l (2.1.12) 

Ker(Q) 



in the sense of infinitesimal calculus of Dieudonne. In what follows, the interpretation of the right 
hand side of the preceding relations will be "a integrated divergence" . In other words, to avoid various 
problems we will use in fact the formal adiabatic limit condition given by ( 2.1.1 1] ) . A more detailed 
discussion on this point can be found in |l7|| . 

By a trivial Lagrangian we mean a Wick expression of the type 

L(x) = d Q N(x) + i-^ I L»(x) (2.1.13) 

with L(x) and L^(x) some Wick polynomials. The first term in the previous formula gives zero by 
factorisation to the physical Hilbert space (according to a previous discussion) and the second one 
gives also zero in the adiabatic limit; this justify the elimination of such expression from the first order 
chronological product T(x). 
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If one completely exploits the condition of gauge invariance in the first order of perturbation theory 
obtaining the generic form of the Yang-Mills interaction of spin-one Bosons up to a trivial Lagrangian. 
We assume the summation convention of the dummy indices a, b, . . . = 1, . . . , r. The result from 
is: 

Theorem 2.1 Let us consider the operator T{x) defined on H-ym as a Lorentz-invariant Wick poly- 
nomial in Aa(x), u a (x), u a (x), $> a ( x ) such that every term has canonical dimension 3 of 4- If it 
verifies the formal adiabatic limit condition then it has, up to a trivial Lagrangian, the following form: 

T YM (x) = f abc [: A^(x)A bl/ (x)d v A^(x) : - : A^(x)u b (x)d^u c (x) :] , 
+f'abc [■ ^a{x)d fl ^> b {x)A^(x) : -m b : <S> a (x)A bfl (x)A£(x) : -m b : $ a (x)u b (x)u c (x) :} 

+f abc : $ a (x)$ 6 (x)$ c (x) : +g abcd : $ a (x)$ 6 (x)$ c (x)$ d (x) : (2.1.14) 

The various constants from the preceding expression are constrained by the following conditions: 

- the expressions f abc are completely antisymmetric 

fabc = —fbac = ~ f acb (2.1.15) 

and verify: 

(m a -m b )f abc = 0, iff m c = 0, Va, b = 1, . . . , r; (2.1.16) 

- the expressions f' ahc are antisymmetric in the indices a and b: 

fabc = -f'bav (2-1.17) 

verify the relation: 

(mf -m b H )f' abc = 0, iff m a = m b = m c = 0, Va,6=l,...,r (2.1.18) 
and are connected to f abc by: 

fab c m c = f' cab m a - f' cba m b , Va, b, c = 1, . . . , r; (2.1.19) 

- the expressions f abc remain undetermined for m a — m b — m c — and for the opposite case are 
given by: 

fabc = -^fabc [{<? ~ «) 2 ~ml + ml] , (2.1.20) 

for m c ytz 0. 

- the expressions g a bcd are non-zero only for m a = m b = m c = rrid = and in this case they are 
completely symmetric. 

Remark 2.2 The presence of indices of type lib and lie is taken into account by requiring that the 
constants from T(x) are null if one of the indices a, b, c takes such values. One can see that this does 
not affect the equations from the statement of the theorem. 

We also have: 
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Corollary 2.3 In the condition of the preceding theorem, one has: 



d Q T(x) = id^{x) (2.1.21) 

where: 

T» = f abc h u a A bu F^ :~: u a u b d»u c A + f' abc (m a : A£$ b u c : + : * B ^*6«c • (2-1.22) 

The expression T(x) from the preceding theorem verifies the unitarity condition 

T(x) f = T(x) 

if and only if the constants f abc , f' abc and f" abc , have real values; it also verifies the causality condition: 

[T(x),T(y)\ = 0, Vx,yGM 4 s.t. (x - yf < 0. 

We close this Subsection with some remarks. 

Remark 2.4 One can see that the necessity of using ghost fields stems from the fact that there seems 
to be impossible to construct the interaction Lagrangian without them. However, from a fundamental 
point of view, one can consider them only as some catalysers [Tft] and hope that one will be able to 
reformulate the whole theory without them. 

Remark 2.5 In the first order analysis one can also use instead of the formal adiabatic limit condition 



(2.1.11) the more physical condition (2.1. It) because no problems connected with the adiabatic limit 
exists in this case. However, as notices in \l3(j, the condition does essentially eliminate the tri-linear 
terms and one looses a lot of the information of the preceding theorem. This is another indication that 
one should work with the formal adiabatic limit condition. 

Remark 2.6 In $ one can find a discussion showing that trivial Lagrangians do not produce effects 
in the higher orders of perturbation theory. 

2.2 Yang-Mills Fields coupled to Matter 

We study here the possibility of coupling Yang- Mills fields to "matter" . We suppose that we are given 
the Hilbert space of "matter" Ti ma tter which is ussually also a Fock space. Then the coupled system 
is described in the tensor product Hilbert space Tym ® Kmatter- One can describe this Fock space 
considering H-ym = ^-ym ® ^matter with the corresponding supercharge operator and forming the 
quotient Ker{Q) / 'lm(Q). We will consider here that the "matter" is formed from Dirac Fermions 
only. 



First, we generalize theorem 2.1 



Theorem 2.7 Let us consider the operator T(x) defined on 'HyM w hich is a Lorentz-invariant Wick 
polynomial in A%(x), u a (x) ,u a (x) ,<3? a (x) and the matter fields such that every term has canonical 
dimension 3 or 4- Then T(x) verifies the formal adiabatic limit condition, if and only if, up to a trivial 
Lagrangian it has the following form: 

T(x) = T YM {x) + A%(x)j aiM (x) + V — *a(aOfytf (a) + Y, *a{x)j a (x) + T mattcr (x) (2.2.1) 

m B ^0 m a =0 
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Here T YM (x) has been defined in theorem 2.1, j afJ and j a are Lorentz covariant currents build only 

from the matter fields with uj(j a n) = 1,2,3 and T ma tter(x) contains only the matter fields. Moreover 
the following conservation law should be valid: 

d^(x)=0, Vm a = 0. (2.2.2) 

The expression for T{x) verifies the unitarity requirement if and only if we have: 

Va = l,...,r, j a (x)t = j a (x), Vm a = (2.2.3) 

and verifies the causality condition if and only if: 

(V)} = 0, (x-y) 2 <0, Va,6= l,...,r, (2.2.4) 

[ja(x),j b (y)}=0, (x-y) 2 <0, Vm a = m b = 0, (2.2.5) 

\j£(x),j b (x)} = 0, (x-y) 2 <0, Vm fe = 0. (2.2.6) 



Proof: 



Beside the terms considered in theorem 2.1 we have to include terms containing explicitly the Dirac 



Fermions. Lorentz covariance and power counting limit these terms to T mat ter{x) and: 

T mat ter(x) = A^(x)j a ^(x) + $a(x)ja(x) (2.2.7) 

with j a ^ {ja) a Lorentz covariant (resp. invariant) operator. Proceeding in the same way as for the 



proof of theorem 2A , we obtain a supplementary restriction, namely: 

maja = dtf£, Va = l,...,r. (2.2.? 

In other words, for m a = we get ([2.2.2 ) and for m a / we get: 

1 



(2.2.9) 



The expression from the statement emerges. The other assertions are straightforward, although 
rather tedious to verify. ■ 

It is clear if the Hilbert space of the matter fields is also a Fock space and the currents are build 
from Wick monomials, then the commutation relations ( |2.2.6j ) are always verified. 

Corollary 2.8 The following formula is true 



d_ 

where T^(x) is obtained by adding to the corresponding expression from the pure Yang-Mills case - see 



d Q T(x) = ig^T^x) (2.2.10) 



(2.1.2i) - the following contribution due to the presence of matter: 

T^ tter {x)^u a {x)f a {x). (2.2.11) 

Now we get in detail the structure of the interaction Lagrangian in the following two propositions. 
We have: 
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Proposition 2.9 Suppose that the Dirac Fermions generating Umatter ar ^ if) a °f niasses Ma > 
0, A = 1, . . . , N . Then the generic form of the currents from the preceding theorem are: 



and 



Ja{x) =■ ^A(x)(t a )AB^B(x) : + : V^XOabT^TsV'.bOzO : 



j a (x) =: 4> A (x)(s a ) AB ipB(x) : + : ip A {x)(s' a ) A BlbipB{x) : 



(2.2.12) 



(2.2.13) 



The causality conditions from theorem \l. % are fulfilled and the hermiticity conditions are equivalent 
with the fact that the complex N x N matrices t a , t' a , s a , a = 1, . . . r are hermitian and s' a , a = 
1, . . . ,r anti-hermitian. 

The contributions with (without) the matrix 75 is called axial (resp. vectorial) current. Let us 
define the mass matrix by: 

Mab = Sa,bM a , VA,B = 1,...,N. (2.2.14) 

Then we have: 

Proposition 2.10 The following mass relations are true: 



-{M,t' a }, Vm a /0, 



[Af,t a ]=0, {M,Q = 0, Vm a = 0. 



(2.2.15) 



(2.2.16) 



In particular, the matrices t a , Vm a = can be exhibited into a block diagonal structure (eventually 
after a relabelling of the Dirac fields) and the masses corresponding to the same block must be equal. 



Proof: It is easy to show that the conservation law (2.2.8) is equivalent to the two relations from 
the statement. ■ 



Corollary 2.11 Let us define 

t e a = t a + et' a , s e a = s a + es' a , Va = 1, . . . , r, e = ±; 



here e = db. Then, the relations (2.2.15) and (2.2.16) are equivalent to: 

sl = —{Mf a -t-*M), Vm a /0, 

Mt e a = t~ e M, Vm a = 
and the hermiticity conditions are equivalent to: 



(2.2.17) 

(2.2.18) 
(2.2.19) 



(*a)*=C «)* = *r, Va = l,...,r, e = ±. 



(2.2.20) 
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3 Perturbation Theory 



3.1 The General Framework 

We give here the basic ideas of a multi-Lagrangian perturbation theory following and fl8|| . One 
can argue that the S-matrix is formal series of operator valued distributions: 

S(g) = 1 + V-r / dx 1 ---dx n T ju ,„ dn (x 1 ,---,x n )g jl (x l )---g jn (x n ), (3.1.1) 
n=1 n\ J R 4 n 

where g = (gj(x)) . =1 p is a multi-valued tempered test function in the Minkowski space M 4 that 
switches the interaction and Tj lr ..j n (x\, ■ ■ ■ ,x n ) are operator-valued distributions acting in the Fock 
space of some collection of free fields with a common dense domain of definition Dq. These operator- 
valued distributions are called chronological products and verify some properties called Bogoliubov 
axioms. We note that there is a canonical association of the point Xi and the index ji. One starts 
from a set of interaction Lagrangians Tj(x), j = 1,...,P and tries to construct the whole series 
n > 2. 

We outline briefly the set of axioms imposed on the chronological products Tj ly „ ; j n ; we do not give 
the explicit formulas because they are well known in the literature and can be found in the references 
quoted above. 

• Symmetry: this axiom describes the behaviour of the chronological products with respect to the 
permutation of the couples (xi,ji). 

• Poincare invariance: this axiom describes the behaviour of the chronological products with 
respect to the action of the Poincare group in the Fock space of the system. Essentially is a 
tensorial covariance condition. 

• Causality: it describes factorization properties of the chronological products for causally sep- 
arated arguments. This seems to be the central axiom of this axiomatic approach; it plays a 
major role in other axiomatic schemes as well. 

• Unitarity: this axiom is considered in the sense of formal series. 

A renormalization theory is the possibility to construct such a 5-matrix starting from the first 
order terms: Tj(x), j = 1, . . . ,P which are linearly independent Wick polynomials called interaction 
Lagrangians which should verify the corresponding axioms expressing the behaviour with respect to 
Poincare transformations, Hermitian conjugation and commutation properties for space-like separated 
arguments. 

The case of a single Lagrangian corresponds to a single coupling constant, that is P = 1 and in 
that case the chronological products will be operators T(X) without any indices. However, it is more 
convenient to consider that the interaction Lagrangian is given by the sum 

T{x) = Y J c j T j {x) (3.1.2) 

with Cj some real constants. In this case, the chronological products of the theory are 

T(X) = ^c n ...c Jn T Ju ... Jn (X). (3.1.3) 
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It can be showed that that one must consider the given interaction Lagrangians Tj(x) to be all 
Wick monomials canonical dimension Uj < 4 (J = 1, . . . , P) acting in the Fock space of the system. 
Because the Fock space is generated by some free relativistic fields acting on the vacuum f2 it is easy 
to see that there always have covariance properties with respect to Poincare transformations. 

If there are non-Hermitian free fields acting in the Fock space, we have in general: 

r j (x)t = r i .(x) (3.1.4) 

where j — > j* is a bijective map of the numbers 1, 2, . . . , P. 

If there are Fermi or ghost fields acting in the Fock space, the causality property is in general: 

T h ( Xl )T h (x 2 ) = {-iyi^T j2 (x 2 )T h {x x ), Vxi ~ x 2 . (3.1.5) 

Here G{ is the number of Fermi and ghost fields factors in the Wick monomial Tj; if o~j is even 
(odd) we call the index j even (resp. odd). One has to keep track of these signs in the symmetry 
axiom for the chronological products. 

It is convenient to let the index j have the value also and we put by definition 

To = 1. (3.1.6) 

Moreover, we define a new sum operation of two indices ji,j 2 = 1,...,P; this summation is 
denoted by + but should not be confused with the ordinary sum. By definition we have: 

T h+j2 (x)=c: T h (x)T j2 (x) : (3.1.7) 

for some positive constant c. We define componentwise the summation for n-tuples J = {ji, . . . , j n }. 
The new summation is non-commutative if Fermi or ghost fields are present. 
We will use the notation 

ojj = J2^j (3.1-8) 

and we call it the canonical dimension of Tj{X). 

Let us denote by u{d) the order of singularity of the numerical distribution d. we use the definition 
from |^] although one can also use the scaling degree introduced by Steinmann (see p0| ). 

Then we add a new axiom, namely the following Wick expansion property of the chronological 
products is valid: 

Tj(X)= Y, e tePO W l(X) (3-1.9) 

K+L=J 

where: (a) tic(X) are numerical distributions (the renormalized Feynman amplitudes); (b) the degree 
of singularity is restricted by the following relation: 

w(tjr) < u K - 4(n - 1); (3.1.10) 

(c) e is the sign coming from permutation of Fermi fields; (d) we have introduced the notation 

Wj(X)=:T jl (x l )---T jn (x n ): (3.1.11) 

Let us notice that from ( 3.1.9| ) we have: 

tj(x) =< n, Tj(x)n > . (3.1.12) 



n 



In particular, these numerical distributions have causal support and are Poincare covariant; trans- 
lation invariance implies that they are in fact distributions in m = 4(|X| — 1) variables. 

The recursive construction assumes that we have the expressions Tj(X) for \X\ < n — 1 verifying 
all the properties and tries to construct them for X = n. The basic object is the commutator function: 

Dji,...,j n ? • • • ; -En—l] %n) — 0^1 > • • • ; ■^n— 1 j ^n) ^ji,---,jn ■ ■ ■ i •^n— 1; •^n) (3.1.13) 

where 

4 1 ,..,in(^>---^n-i;x n )^^ , Xi)X26pari{x) {-l)\ x >\ Tjl (Xi)f*{X 2 ) (3.1.14) 

and 

Z^...^,...,^;^ (3.1.15) 

the sums ^ run over the partitions verifying X2 7^ 0, a; n £ X\. 

The commutator function can be proved to be Poincare covariant and to have causal support i.e. 
supp(Dj lj ^,j n (xi, . . . , x n -i; x n )) C T + (x n ) U T~(x n ) where we use standard notations: 

T ± (x n ) = {(x 1 ,...,x n )e(M. 4 ) n \x l -x n €V ± , Vi = l,...,n-1}. (3.1.16) 

Moreover, a formula similar to ( 3.1,9| ) is true: 

Dj(X)= e d K{X) W L (X) (3.1.17) 

K+L=J 

where dx(X) are numerical distributions; in analogy to ( |3.1.12| ) we have: 

dj(x) =< n, Dj(x)n > . (3.1.18) 

It follows that the numerical distributions dj{X) have causal support i.e supp(dj{X)) C T + (x n ) U 
T~(x n ) and are SL(2, C)-invariant. Moreover, their degree of singularity is restricted by 

uj(d K ) < uj k -4(n- 1); (3.1.19) 

this is the content of the power counting theorem. One knows that there exists a causal splitting 

dj = aj - rj, supp(aj) C T + (x n ), supp(rj) C T~(x n ) (3.1.20) 

which is also SL(2, C)-invariant and such that the order of the singularity is preserved. So, there 
exists a SL(2, C)-covariant causal splitting: 

Dj(X) = Aj(X) - Rj(X), \X\=n (3.1.21) 

with supp(Aj 1> ___j n (x 1 , . . . ,x n -i;x n )) C T + (x n ) and supp(R jl) ...j n (xi, . . . , x n _i; x n )) cT~(x n ). 
Let us define 

Tj(X) = Aj(X) - A'j(X) = Rj(X) - R'j(X). (3.1.22) 

Then these expressions satisfy the SL(2, (C)-covariance, and causality axioms. One can also fix uni- 
tarity and symmetry. 

We end this Subsection with an important remark. Let us consider some general Wick polynomials 

A i (x) = J2 c ij T i( x )i 2 = 1,2,... (3.1.23) 
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Then we can define the chronological products: 

T(Ai( Xl ),--- ,A n (x n )) = ^2c hjl ---Ci njn Tj lt _j n (x\, ■ ■ ■ ,x n ). (3.1.24) 

J 

One can find in [ 10 1 a system of axioms for the expressions T(Ai(xi), • • • , A n (x n )) which is equiv- 
alent to the Bogoliubov set of axioms. 

3.2 Ward Identities 



As we have said in the Subsection 2.1 the problem is to construct the whole series T(X) such that 
one has the gauge invariance condition in all orders of the perturbation theory at the same time with 
the other Bogoliubov axioms. 

In general we have something more general than relation ( |3.1.2| ): 

T{x) = Y J C 3 T j (x) T"(z) = £c£ZK*) (3.2.1) 

with cj, Cj some real constants; then we will have something more general than fl3.1,3|) : 

T(X) n j h jn (X), Tf{X) Cj n Tn,..., jn ( X ) (3-2-2) 

In particular, the following conventions hold: 

T(0) = 1, 2f(0) = 0, If (X) = 0, for x t X. (3.2.3) 
Then the gauge invariance condition ( [2.1.11] ) can be written more compactly as follows: 

d Q T{X)=iJ2^pT?{X). (3.2.4) 
We suppose that these relations are true up to order |X| < n — 1 and investigate the possible 



obstructions in order n. The procedure used in flilfl , [jl^ l and Hxql , [17] amounts to the following. Let 
us define the operator distributions D(X) and D^{X) in analogy to the relations ( |3.2.2| ). Then it can 
be proved that we have: 

n Pi 

d Q D(X) = i g^ D i( X )^ \ X \ = n - ( 3 -2-5) 
i=i %l 

We can express this condition in terms of numerical distributions. According to the relation ( |3.1.S| ) 
and Wick theorem we must have Wick expansions for the two expressions appearing in the preceding 
equation: 

D{X) = Y J dj{X)Wj{X) : D?(X) = Y,d?.j(X)Wj(X) (3.2.6) 
J J 
The numerical distributions appearing in these relations have the following properties: they are 
Poincare covariant, they have causal support and the order of singularity can be restricted according 
to the power counting formula: 

uj{d,j) + loj < 4, u(df.j) +ujj<A (3.2.7) 
according to the power counting theorem. 
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One can rewrite ( |3.2.6| ) as follows: 

D(X) =Y^d i (X)W i (X), D»(X) =Y,di{X)W^X) + ^d? (X)W l;i (X) (3.2.8) 

i i % 

where di and df can be taken linear independent over the vector space of smooth functions with 
polynomial bounded increase at infinity Om- The index i takes a finite number of values and the 
expressions Wi(X), Wi-i(X), W^X) are Wick polynomials. 

Using the linear independence one obtains from ( 3.2. 5| ) a set of identities among Wick polynomials 
of the type 

dqWi = ■■■ (3.2.9) 

where the left hand side can be computed as follows. First one makes the derivation operations in the 
right hand side of (|3.2.5|) . It is quite possible that relations of the type 

^d».(X)=Y,c j d j {X) (3-2.10) 
1 j 

are valid for some numbers c,-. Then one has to rearrange the expression in the right hand side of 
(|]2j) and the right hand side of ( ^9|) emerges as the coefficient of di(X). 

Identities of the type ( 3.2.10[ ) are called Ward-Takahashi (or Slavnov- Taylor identities). In [12 
these relations are called the C-g identities. They have been extensively studied in || . In lower orders 
of perturbation theory one can check them by explicit computation. 

One now can interpret the renormalization theory as a distribution-splitting preserving of the Ward 
identities. Suppose that one can find a causal splitting di = df dv — d\ et of the set of causal distributions 
di(X) such that we preserve Poincare covariance, the order of singularity and the identities ( |3.2,10| ) 
i.e. we also have: 

_^_ (d M )«*,(«*)(*) = ^ Cj d a j dv{ret) (X) (3.2.11) 
1 j 

Then we define the expressions A(X) and A^(X) by making into the formulae ( |3.2.8| ) the substi- 
tutions d — > d adv . If we use now the relations fl3.2.9| ) we easily obtain 

n f> 

d Q A(X)=iJ2g^Af(X), \X\=n. (3.2.12) 
i=i Xl 

The similar property for the chronological products of order n easily follows. So, the obstructions 
to the gauge invariance in order n can appear in the process of causally splitting the relations Q3.2.10| ) 
i.e. we might have instead of (3.2.11): 

±. {d ^r d ^{X) -Y.Cjdf^iX) =p(X) (3.2.13) 

1 3 

where the expression in the right hand side p{X) - called anomaly - must have the form 

p(X) = p(d)5(X); (3.2.14) 
here p(d) is a Lorentz covariant polynomial in the partial derivative operators and 

5(X) = 5(x 1 - x n ) ■ ■ ■ 5(x n -i - x n ). (3.2.15) 
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Also, if the distribution appearing in ( 3.2. 10| ) have some global symmetry property (symmetry with 



respect to some global group of symmetries, (anti) symmetry with respect to some indices, etc.) one 
can usually perform the distribution splitting such that these properties are preserved also. Moreover, 
we have a limitation on the degree of the polynomial p(d): 

deg{p) < u (3.2.16) 



where u is the order of singularity of the left hand side of ( 3.2.13| ). It easily follows a case when there 
are no anomalies, namely when: uj(df) < 
of a anomaly can be simplified by redefini 
adding expressions of the type p(d)5(X). 



are no anomalies, namely when: ui{d^) < —2, V//. Let us note in closing this Section that the form 



of a anomaly can be simplified by redefinitions of the distributions aj and a? ■; we have the freedom of 
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4 Second Order Perturbation Theory 
4.1 Yang-Mills coupled to Matter 

We follow [|l7]] where the pure Yang-Mills case was studied emphasizing the possible appearance of 
anomalies in a more explicit way. 



Theorem 4.1 Suppose that the distribution T(x,y) verifies ( 3.1.10J . The it verifies the formal adia- 
batic limit condition if and only if the following identities are verified: 



fabcfdec fbdcfaec "t~ fdacfbec — 0, O, b, d, e — 1, . . 
fdcafceb ~ fdcbfcea = ~~ fabcfdec a,b, d, C = 1, . . . 
fcabfcde + fcdbfcae + fcebfcda = 0; ^ m a = 171b = m d 
Sbcdeffcbadcdef = 0, a, b, d, e, f = 1, . . . , r\ 

[*o> *&] = ifabctl, e = ±, a, 6=1, ...,r; 
*o S ft~ ~~ S b~*a = ifbca s t ■> a, 6 = 1, . . . , r. 



. r; 



0; 



(4.1.1) 
(4.1.2) 
(4.1.3) 
(4.1.4) 
(4.1.5) 
(4.1.6) 



Here S... is the symmetrization operator in the indices which are explicitly exhibited. 
Proof: (i) According to the ideas from Subsection |3.2| , we compute the commutator 

D( Xl ,x 2 ) = [T{xi),T(x 2 )] 



(4.1.7) 



using Wick theorem and identify a set of linearly independent distributions d{ as in (3.2.8); these 
are distributions in one variable £ = x\ — x 2 due to translation invariance. Direct inspection of the 
expressions ( p. 1.14 ) and ( 2.2.1] ) produces the a list of such distributions A with causal support. Using 
Feynman graph terminology we have distributions associated to tree and one, two and three loops 
graphs. All these distributions can be written as sum of the positive (negative) frequency parts: 



(a) From tree graphs: 



a = a(+) + aH 



n f) n 8 8 D 

S M (x) = (i>yd + M)D M (x), 



(4.1. 



(4.1.9) 



where D m is the Pauli-Villars commutator distribution of causal support corresponding to mass m 
(see [|l7|] for the definition) and Sm is the similar distribution for a Dirac field of mass M. 
(b) From onedoop graphs we get new distributions with causal support. 



D (±) 

"mi ,m2 



D (±) 
mi,m 2 ;p 



±D$(x)dW(x), 
±D^d p D% - (1 - 2). 
d„D, 



•^mi,m2;pff 



8 nWfl nW - nWfl /j n(±) 

u P ly m 1 u cr 1J m2 "mi u p u o LJ m 2 



p'-'m\ ,ni2 j 

+ (1^2) 
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(-), 



?(+)/ 



P 



(±) 



M 1 ,M 2 ;p 



(x)^±Tr s[- ) i (Txh P S^>(± 



(+), 



3 (±) 

M\,Mi,pa 



(x)^±Tr lp S\->(Tx) 7a S i + ) 2 (±x) 



J m,M 



+ n(±) o(±) 



(4.1.10) 



We note that in the definition of Dm},m 2 ;p we have taken the antisymmetric part in the masses 
because the symmetric part has been considered separately :it is the third distribution from the list, 
(c) From two- loops graphs: 



" mi, ma, wis 



D (±) D (±) D (±) 
"mi m 2 7713 ; 



<9 2 P, 



"mi ,ma ;m3 



P. 



mi ,m2 .m3 
u p. LJ ra\ u ^m2 "7773 ) 

= +d( ± )p (±) 

m;Mi,M 2 ~ r M\,M?. > 

= ±n(±)p( ± ) 



P. 



(±) 



m;Mi,M 2 ;pcr L "m - 1 Mi,M 2 ;pcr 

(d) From three-loops graphs: 

L» (±) fx) = iljW D {±) D (±) 

"mi ,m 2 ,7773,7774 V'*',' — -"-"mi "m 2 7773 "m4 



(4.1.11) 



(4.1.12) 



The distributions P '■ are obtained from contractions of two vectorial currents. Let us note that one 
also obtains distributions of the type Q \ from contractions of two axial currents. These distributions 
can be obtained directly from the corresponding distributions P "' by conveniently inserting two 75 
factors. However, the distributions of the type Q " can be expressed in terms of P'" if one uses the 
identity 

75-Sl? 75 = (4.1.13) 

The distributions following from contractions of axial and a vectorial current are null because the 
traces so obtained are null. Next, we note that in the other commutators 



Dl( Xl ,X2) = [T»(xi),T(x 2 )}, D2( Xl ,x 2 ) = [T( Xl ),T»(x 2 )} = 
the distributions g^xd^ from ( |3.2.8|) can be of the following type: 

dpDmi ^/pSpj, -Dmi,m 2 ;/75 ^mi,m2;/ii/) PM\,M2\pi 

and the distributions of the type di can be of the type: 



-Df(x 2 ,xi) 



Mi,M2\pv 



D r 



8 D 



D 



mi,m 2 1 



D 



m\m,2\p- 



(4.1.14) 
(4.1.15) 
(4.1.16) 



Here the various parameters m, M, . . . are the masses appearing in the theory. If we consider 
distinct combinations of masses and indices we obtain a linear independent set of distributions. 

Let us also give for further use the orders of singularity of the distributions listed above. We have: 

w(An) = —2, u(D mijTn2 ) = 0, Lo(D mi ^ m2 - p ) = —1, Lo(D mi ^ m2 - pa ) = 2, w(Pm 1 ,m 2 ) = 2, 

w(PMi,M 2 ;p) = 1, w(Pm 1 ,M 2 ;p ( t) = 2, Uj(P m -M 1 ,M 2 ) = 4, w(P m; Af 1 ,M 2 ;pcr) = 4, 
m,Mj — 1] ^ (-^mi ,7712 ,m3 ) — 2, CL>(-D miim2 ;7773) — 4, ^(-^mi, 777^7773,7774) — 

4. (4.1.17) 



17 



Some of these orders of singularity are in fact lower than naive power counting suggests. 
All these distributions have causal support so we have causal decompositions: 



A = A adv - A ret (4.1.18) 

We have assumed that the causal splitting is preserving Lorentz covariance and the order of 
singularity. If the order of singularity is less then this causal decomposition is a unique (see the end 
of the preceding Section). This is the case for the distributions D m , Sm and D mim2 - p . 

(ii) Now we consider the Ward identities ( |3.2.10| ). By direct inspection one finds out that they are: 

(d 2 + m 2 )D m = 0, (4.1.19) 

(*7 • a - m)s m = s M (ir a -m) = o, (4.1.20) 

ml)D mi>m2 , (4.1.21) 

■m{)D miim2;v (4.1.22) 

-M 2 )P Ml , M 2 (4-1-23) 

- M 2 )P Ml ,Mw (4.1.24) 

Now we analyze possible anomalies resulting after the causal splitting procedure. It is well known 
that the first two relations ( 4.1.19Q and ( |4.1.20| ) are indeed producing anomalies: for the (unique) 





= (ml 




= ( m 2 


9' 1 Pm 1 ,M2;^ 


= i(Mi 




= i(Mi 



causal splitting considered above one gets: 

(8 2 + m 2 )D^'< ret '> =5, (4.1.25) 

(i 7 • a - M)slf {ret) = sfj" {ret) (i T a -m) = is. (4.1.26) 

One can prove more than that: even if we modify these splitting with arbitrary local polynomial 
terms the anomalies do not disappear. 



Next we consider ( 4.1.22 ); inspecting the orders of we can have the following generic form of the 
anomaly: 

Pu(d) = c 1 8 u + c 3 d l/ d 2 . (4.1.27) 
We can eliminate this anomaly if we make the redefinition: 

Dt,mr,^ - £Cm 2 ;^ + (ci^ + ^8 V )S. (4.1.28) 



The case ([4.1.24 ) can be treated in a similar way and the anomaly is also eliminated. The Ward 



identity ( 4.1.2l| ) is non-trivial only for m\ ^ We have already noticed that there exists a unique 



causal decomposition preserving Lorentz covariance and the order of singularity of D mij7nr ^\ then we 
can define: 

^mi,tn 2 = 2 2^^ m l> m 2;M (4.1.29) 

' m-2 — m| 



and the relation ( 4.1.21| ) is preserved; moreover the order of singularity is preserved: ^{Dmim 2 ' 



The Ward identity ( |4.1.23| ) is non-trivial only for Mi ^ M 2 and it has the generic form 

p(8) = c + c 2 d 2 . (4.1.30) 
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If we make the redefinitions: 

Pmi,m 2 ;h ~* Pmi,m 2 -,h + c 2dfj,5, Pmi,m 2 ~* Plh,M2 + ^ m i — M 2 ^ (4.1.31) 

the anomaly is eliminated. 

It is interesting to summarize the preceding argument by saying that the anomalies are produced 
only by the distributions associated to tree graphs. 

(iii) It follows that we can describe the structure of the terms from D p (xi, x 2 ) which can produce 
anomalies. It is sufficient to consider I = 1 and notice that the other piece doubles the value of the 
anomaly (because of obvious symmetry properties). We have: 

d d 2 
Df(x 1 ,x 2 ) = -^—D mc (xi - x 2 )T c (x 1 ,x 2 ) + dx p 1 Dm ^ Xl ~ X2 ) T c( Xl ' X2 ) 

d , d 2 

+-^—D m * c (x 1 - x 2 )T c (xi,x 2 ) + g D m * (x 1 - x 2 )T'/{xi 1 x 2 ) 

8 8 

+ Y, ■ U^\x l )rS A {x l - x 2 )V { «\x 2 ) + ■ T i / ) (x 1 )S A (x 1 - x 2 )^W ( A a) (x 2 ) : + • • • (4.1.32) 

a=l a=l 

where by • ■ ■ we mean the contributions which do not produce anomalies because of the argument of 
(ii). We have the following explicit expressions: 

T c (xi,x 2 ) = T^ M (xi,x 2 ) + f abc : u a (xx)A p (xi)j cp (x 2 ) : 

T' c (xi,X 2 ) = T c ' yM (xi,X 2 ) - f'cab : $a{xi)u b {xi)j c {x 2 ) \ 

T p (x x ,x 2 ) = T c yM ^(x!,x 2 ), T'J>{x x ,x 2 ) = T? M >P( X1:X2 ) (4.1.33) 

where the expressions Tj M (x\,x 2 ), T^ m (x\, x 2 ) and tJ M ' p {x\, x 2 ), Tc YM,p (x\,x 2 ) corresponds 
to the pure Yang-Mills and can be found in p7|. Also 



Uf{x) = Uf{x) = Uf\x) = U A 7 \x) = (t a ) BA u a (x)ii B (x), 
Uf{x) = Uf{x) = Xjf{x) = U A 8 \x) = -(t' a ) BA u a (xj$ B (x)j 5 , 
^a\ x ) = Va ] ( x ) = {t h ) ADl ^ D {x)A p {x), V [ A 2 \x) = Vf(x) = -(t' b ) ADlpT My)A p b (x), 
Vf(x) = Vf(x) = (s b ) AD ^ D (x)^ b {x), Vf\x) = V (7) (x) = {s' b ) ADl ^ D (x)^ b (x). (4.1.34) 



and 



W%\x) = Wf\x) = Wf{x) = W A 7) (x) = (t a ) AB u a (x)ij B (x), 
W A 2 \x) = wf{x) = wf(x) = Wf{x) = {Q BA u a {x)^ B {x\ 



A K^J — yy A W ~~ rr A ~ A 



TX>(x) = T™(x) = -{t b ) CA ^ c (x) lp A p {x), TX>(x) = T A *>(x) = -(t' b ) CA ^hp^A p (x), 
Tf{x) = Tf{x) = -(s b ) CA ^ c (^b(x), Tf\x) = T A 7 \x) = -{s' h ) CA ^c^)l^b{x). (4.1.35) 

The expression of the anomaly can be obtained in the generic form: 

A( Xl ,x 2 ) = 6(xi - x 2 )A{xi) (4.1.36) 
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where: 



C 



T c (xi,xi) + I^(xi,aci) 



_5 



p 2£ ) (zi,a:i) 



p^ P ) (xi,Xl) 



+ 



(4.1.37) 



So, the expression of the anomaly A(x) gets an extra term because of the presence of the Dirac 
Fermions: 

A(x) = A YM {x) + i : u a (x)A^(x)^ A (x)j p ([t a ,t b ] + [t' a ,t' b ] - if a bct c )AB^B(x) 
+i : u a (x)A^(x)i/j A (x)-f p -f 5 ([t a ,t' b \ + [t' a ,t b ] - ifabct'c)ABipB(x) 
+i : u a (x)^ b (x)ip A (x)([t a , s b ] - {t' a , s' b } + if' cba s c ) A B^B(x) 
+i : u a (x)^ b (x)ip A (x)^ 5 ([t a , s' b ] - {t' a , s b } + if' cba s' c )AB^B{x) : (4.1.38) 

(iv) We proceed now as in jL7]]. First we equate the expression A{x) to a coboundary oIqL(x). 
We get all the relations from [|l7j (and this explains the first four relations from the statement). 
Moreover we obtain for all a, b = 1, . . . , r: 



[ta,S b ] ~ {t' a ,s' b } = -if' c 



cba a ci 



[ta, S b ] — {t' a , S b } — —if' cba S c . 



(4.1.39) 



which are equivalent to the last two relations from the statement. 

(v) From the preceding computations we can obtain the explicit expression for the coboundary 
L(x): it coincides with the expression obtained for the pure Yang-Mills case: 



L(x) = L 



YM, 



x) = -^fcabfcde ■ A av (x)A bv (x)A^(x)A V e {x) 
-fcdafceb ■ A av {x)Al{x)<S> d (x)<S> e {x) 
~ Yl 9abcd ■ ®a{x)$ b (x)®d{x)<S> e (x) 



m b ^0 



where: 



9abcd 



1 



'Sabdefcabfcde- 



Let us also define: 



2m b 

L^x) = J2[T c ^^)+T^^)] 

c 

Again it coincides with the expression from the pure Yang-Mills case: 



(4.1.40) 

(4.1.41) 
(4.1.42) 



L"(x) = L™*{x) = -f cab f cde : u a (x)A bu (x)A» d (x)A£(x) : -f cab f cde ■ ® a {x)u b (x)<5> d (x)A»{x) : 

(4.1.43) 

We consider now a canonical causal splitting A c (xi,x 2 ) and A?'^(xx, x 2 ) given by the expressions 
which are obtained from the corresponding commutators if we make the substitutions: A — > /\ adv . 
This indeed gives a causal splitting of D{x\,x 2 ) and D^(xi,x 2 ) respectively. However the identity 
d3T27l2|) is not fulfilled. If we define now the new causal splitting 

A(xi,x 2 ) = A c ( Xl , x 2 ) + S(x! - x 2 )L(xt), Af( Xl ,x 2 ) = Af 1 (x 1 ,x 2 ) + 5(xt - x 2 )I/ l (x 1 ) (4.1.44) 
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then one can see that ( |3.2.12j ) becomes true. Moreover, in this way one can obtain in the usual way 
the expression of the chronological products T(x\,X2) and T^(xi, X2) such that we have ( 3.2.4 ) and 
all other properties, in particular symmetry. ■ 

Remark 4.2 If we do not require that ( 3.1.1C ) is fulfilled, the relations ( \4-l-4 ) an d (U-l-Q ) acquire a 
weaker form. 

The group-theoretical informations contained in this theorem are: 

(a) The expressions f abc are the structure constants of a Lie algebra g. 

(b) The structure constants f abc corresponding to m a = m b = m c = generate a Lie subalgebra 
00 C 0. 

(c) The r x r (antisymmetric) matrices T a , a = 1, . . . , r defined according to 



{T a )bc = —f'b, 



Va, 6, c = 1, 



(4.1.45) 



are an r-dimensional representation of the Lie algebra g. 

The representation T a exhibited in the statement of the theorem is nothing else but the represen- 
tation of the gauge algebra g into which the Higgs fields live. 

(d) The relation (4.1.5) tells that the matrices t € a are representations of the Lie algebra g and 
relation ( 4.1.6| ) shows that the matrices s e a are some tensor operators with respect to the couple of 
representations t\ of the Lie algebra g. 

So, we propose the following strategy of analyzing the generalization of the standard model de- 
scribed in this paper: first one should find out restrictions on the Lie algebra g from the relation 
(gXg ), then one takes a couple of representations t e a of this Lie algebra and afterwards one determines 
the matrices from the relation ( [4.1.6 ) using ideas from the proof of Wigner-Eckart theorem. We 
mention that if one tries to substitute the formula (2.2.18) into the formula (4.1.6), as it is done in Q, 
then we end up with some very complicated trilinear relations which are extremely difficult to analyze 
in the general case. 

Next, we have a generalization of Proposition 3.9 from [17]. By definition the Feynman propagator 
and the Feynman antipropagator are: 



Then we have: 



ah 



. ret 



+ A 



(+) 



.AF 



, adv 



, ret 



AH. 



(4.1.46) 



Proposition 4.3 Suppose that that there is no contribution Ti j?natter in the first order chronological 
product. Then, we have 



T c (x,y) 



^YM,c, 



-fabcD^ c (x - y)[: A av (x)F b yp (x)j cp (y) : - : u a (x)d p u b (x)j£(y) : +(x <-> y) 



(x-y)[: AP a {x)A»{x)j cp {y) : -{x 



f —n F 

c dx^ 

-f ahc D F mc (x - y)[: ^ a {x)d^ b {x)f c {y) 
-f abc D^(x - y)[. d^ a (x)A^(x)j c (y) 

-f abc JL D ^{x-y)[. * a {x)A*{x)Mv) : -(.<■ 



+{x *- 


^ y)] 


— (X<r- 


y)] 


— (x <- 


■* y)\ 


+(x *■ 


^ y)\ 


— (X <r 


■* y)] 


+(x *■ 


■* y)] 
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+hZD^ c (x - y)[: A a ^x)A^(x)j c (y) : +(x <- y)\ 
+ h< cl D m* c ( x ~ y)i : u a {x)u b (x)j c {y) : +(x y)] 
+3h^ c D^(x - y)[: * a (x)* b (x)j c (y) : +(x <- y)] 
+Ag ahcd D^ n ,(x - y)[. ® a (x)<f> b (x)<f> c (x)j c (y) : +(x y)] 
+ : A%{x)A p b (y) : {[(t a ) AC {t b ) C B ■ ^ A { x )lA c {x - yhpMv) ■ 

+ (t'a)Ac{t' b )cB : ^a{ x )1^ S M c ( X ~ y)lpl^B{y) ■ 

+{t a ) A c{tb)cB ■ ^a( x )iij,Sm c ( x - y)i P i^B{y) ■ 

+ (t' a )Ac{tb)cB ■ i^Aixh^S^ix - y)j p ^B(y) ■ -{a <-> b,fl <-> p,x <-> y)] 

+ (^a)AB(^)iM^M A M sW >(2 ; ~ I/) + (O^B (4)bA<2m a M s;w >(z ~ 2/)} 
+ : $ a (x)<S> b (y) : {[(s a )Ac(sb)cB : ^aO*0<Sm c (x - y)Y>fi(y) : 
+(*o)ac(4)cb : V'A(^)755 , M c (a; - y)75^s(y) ■ 

+ (s a )Ac(s' b )cB ■ ^a( x )Sm c (x ~ y)75^fi(y) : 

+ (- s a)Ac(sfe)ciJ : V'a(^)75'S'm c ( :c - y)^B{y) ■ ~{a <-> b,x <-> y)] 

+ (s a )AM s &W P M A ,M s O* ; - y) + (4) AB (4)baQm a ,M b ~ 2/)} 

+ : ^(x)^(y) : {[(t a )Ac(s&)cB : V^h/^M c (z - y)Y> B (y) : 
-(s^AcCOcs : ^A(y)SM c (y ~ x)i^b(x) : 
+(*a)Ac(4)cB : ^ a 0*07^75 ^(a; - 2/)75^s(y) : 
-(4)ac(Ocb : V'A(y)755 , M c (y - x^^six) : 

+(* a )Ac(4)cB : V ; a(^)7^5'm c (^ - y)ibi>B{y) ■ 

-(s'b)Ac(t a )cB ■ ^A{y)lbSM c {y - x)j^ B {x) : 
+(ia)Ac(s 6 )cB : ^aO)7^Sm c (x - y)j 5 ^B(y) : 
-(s&)ac(Ocb : ^A(y)755 , M c (y - x^^six) : 

+ (ta)AB(s b )BAPM A ,M B -A X ~ V) + (O As(4) BaQm a ,M B -A X ~ 2/)] ~ !/]} 

- D m a ( x ~y) ■ jan(x)j£(y) : 

-(to)i4c(*o)cB[: V^h^^McO - y)^^ B {y) ■ +(x y)\ 

-(t' a )Ac(t' a ) C B[- *l>A( x )7^5Zm a ,M c ( x ~ vWlb^B (y) ■ H x V)] 
-{ta)Ac(t' a )cB[- ^ A{ x )l^m a ,M c ( x ~ y)l^l^B{y) ■ +{ x V)\ 
-{t'a)Ac{ta)cB[- i> a( x )1 nl^rn a ,Mc ( X ~ f)V^B(j/) = H x V)] 
-9^[(ta)AB(t a ) B AP^ a] M A ,M B ;^( x ~ V) + (QaB^QbaQ^-MaMb^ ~ V)\ 

+ D m*( x ~ y) : ja(x)ja(y) ■ 
+ {s a )Ac(s a )cB[- i>A( x )^m a ,M C ( X ~ P^fo) ■ + ( x V)] 
+ (s' a )Ac(s' a ) C B[- V> ' A( x )l^m a ,Mc( X ~ y)75^fi(y) : + (& y)\ 
+ (s a )Ac(s' a )cB[- tpA( x )^m a ,M c ( x ~ y)lb^B{y) ■ +(s ^ y)] 

+(Oac(s )cb[: ?a(^)7 5 S^ >Mc (x - y)V's(y) : +(z <- y)] (4.1.47) 
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Here h< abc = Wbca m b + facb m a), = f'abc m b- 

Let us note that the expressions (p. 2.12| ) and (|2.2.13|) for the currents can be also written as follows: 



f a {x) =: ^\{x){tt)ABl^t{x) : + : ^{x){t') ABl ^{x) : (4.1.48) 

and 

j a {x) =: ^{x){st) A B^t{x) : + : ^(x){s~) AB ^(x) : (4.1.49) 

where we have defined for 

= e = ± (4.1.50) 

and the components corresponding to the signs + (resp. — ) are called chiral components of the 
currents. 

4.2 The Conservation of the BRST Current 

The expression 

3brst( x ) = (d-A a + m a <S> a ) T u a (4.2.1) 
is called the BRST current. One can verify easily the conservation of the BRST current: 

^Jbrst = 0- (4-2-2) 

Formally, the BRST charge is given by 

Q = f d 3 xj° BRST (x). (4.2.3) 

We want to investigate the conservation of this current in higher orders of perturbation theory. 
We present here the analysis in the second order. First we have: 

Proposition 4.4 The following relation is verified: 

[jBRST( X ±)> T ( X 2)] = D m a (xi ~ X 2 )A P {x\ , X 2 ) + <9 M An a {x\ ~ X 2 )B a {x\ , X 2 ) 

+d p D ma (x 1 -x 2 )A^(x 1 ,x 2 ) + d"d p D ma (x 1 -x 2 )B p (x 1 ,x 2 ) (4.2.4) 



where 



B a ( Xl ,x 2 ) = hf} c : d u A u a ( Xl )$ b (y)u c (x 2 ) : +m a hf} c : <S> a (xi)$ b (y)u c (x 2 ) : 
- m af' a bc '■ u a {x 1 )dp^ h {x 2 )A p c {x 2 ) : -m a h^ c ■ u a {xi)A hp (x 2 )A p c (x 2 ) : 

(2) " 

' m a h abc '■ u a(xi)u b (x 2 )u c (x 2 ) : -3m a f bca : u a (x)$ b (y)® c (y) : -m a : u a (xi)j b (x 2 ) : (4.2.5) 



and 



B p (x 1 ,x 2 ) = f bca : d u Aa(xi)A^(x2)Uc(x 2 ) : -m a f bca : ^ a (x 1 )A p (x 2 )u c (x 2 ) : 
+fabc ■ u a (x 1 )A bu (x 2 )F^ p (x 2 ) : -fabc ■ u a {xi)u b (x 2 )d p u c (x 2 ) : 
+fi a : u a ( Xl )<S> b (x 2 )d p <S> c (x 2 ) : +m c f' cba : u a ( Xl )^ b (x 2 )A p (x 2 ) : + : u a { Xl )j p (x 2 ) : (4.2.6) 

The expressions for h^J b and h^ c have been given in the preceding Proposition. 
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The computations are long but straightforward. Applying the procedures of the preceding Sub- 
section we obtain from here: 

Proposition 4.5 The following expression 

T c (j^ RST ( Xl ),T(x 2 )) = D^ a ( Xl - X2)A^ Xl ,x 2 ) + d^D^ a ( Xl - x 2 )B a ( Xl ,x 2 ) 

+&>D$ ia {x x -x 2 )A^( Xl ,x 2 ) +d^d p D^ a (xi -x 2 )BP( Xl ,x 2 ) (4.2.7) 

is valid for the canonical chronological product. 

We have the following result which can be interpreted as a conservation of the BRST current in 
the second order of perturbation theory: 

Theorem 4.6 There exists a finite renormalization such that one has the following conservation law: 

jrT(f BRST ^i),T(x 2 )) = i5( Xl - x 2 )^T^x 1 ) (4.2.8) 



Proof: We start from the obvious relation 





dxt 



[(iW(*i)> T (*2)]=0 (4.2.9) 



and perform the canonical causal splitting using the expression of the commutator derived above. If 
we proceed in analogy to the derivation of the consistency conditions for the second order chronological 
products we obtain: 

iT c {f BR sA x i)^2)) = -ijrr [6{x! - x 2 )N»( Xl )] - i5{ Xl - x 2 )A{ Xl ) (4.2.10) 



where: 



and 



(4.2.11) 



iV"(xi) = ^B^(xi,xi). (4.2.12) 

a 

Now it is a matter of computation to prove that we have: A = —d^T^. If we perform the finite 
renormalization T(j RRST (xi),T(x 2 )) = T c (j BRST (xi),T(x 2 )) — i5(x\ — x 2 )N^(xi) then we obtain the 
conclusion from the statement. ■ 

We remark that one can obtain from the preceding result the gauge invariance condition ( [2.1.11 ) 
for n = 1 using the method of Appendix B of [10|. 
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5 Third-Order Gauge Invariance 
5.1 The Derivation of the Anomaly 

In this Section we will analyze the possible obstructions to factorization of the S"-matrix in the third 
order of the perturbation theory. In principle, there is no difference with respect to the preceding 
Section. Nevertheless, the details of distribution splitting are considerably more complicated and the 
same is true for the whole combinatorial argument. 

First we first give a standard regularization procedure of the distributions appearing in the lists 
( 4.1.S| )-( [4.1.12 ). We we choose uiq > different from all masses of the model and write the Pauli- Villars 
distribution for any mass as follows: 

An = D mo + D re °; (5.1.1) 
one can check that the order of singularity of D reg is 

uj(D re9 ) = -4. (5.1.2) 

The decomposition ( pM.l|) induces a similar decomposition for all distributions in the lists ( |4.1.9[ )- 



(4.1.12): we have a sum of two pieces: 

A = A + A reg (5.1.3) 

where 

uj(A°) =w(A), io(A re9 ) = u(A) - 2 (5.1.4) 

and the support properties of A in the momentum space are more convenient. We have: 

A°^(p)~6(±p )f(p 2 ) (5.1.5) 

with supp(f) C {p 2 > A 2 } for some parameter with mass significance A > 0; (for the distributions 
A^(p) we can have in principle A = 0.) 

The main result is contained in the following theorem: 



Theorem 5.1 Suppose that the distribution T(xi,x 2 ,x 3 ) verifies the condition ( 3.1. It ). Then it 
verify the formal adiabatic limit condition if and only if, beside the conditions from the statement of 
theorem \2. | , we also have the following set of supplementary conditions: 

Tr (tt{tt,t+}) = Tr (tj{^,^}) , (5.1.6) 

f'abcd'bfgh + ffbcdbagh + f'gbcd'bafh + fhbcd'bafg = 0- (5.1.7) 



Proof: 

(i) As before, we will investigate the third order commutators 

D(xi, x 2 ; x 3 ) = [T(x 3 ), T(xi, x 2 )} - [T(xi, x 3 ), T(x 2 )] - [T(x 2 , x 3 ), T(xi)] (5.1.8) 



and 



D%(xi,x 2 ;x 3 ) = \T{x 3 ),Ti{xi,x 2 )} - [Tf (xi,x 3 ),T[x 2 )] - [T(x 2 ,x 3 ),T f ^(x 1 )], 
D^(xi,x 2 ;x 3 ) = [T(x 3 ),T2(x 1 ,x 2 )] - [T(xi, x 3 ), T± {x 2 )\ - [T^(x 2 ,x 3 ),T(x 1 )], 
D$( Xl ,x 2 ;x 3 ) = [Tf(x 3 ),T(xi,x 2 )] - [T»{ Xl , x 3 ), T(x 2 )) - [T^(x 2 ,x 3 ),T{x 1 )}. (5.1.9) 
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All these operator-valued distributions have the causal support property, 
(ii) We convene to denote by generically by 



A^ +) (xi - x 2 ) = \\ < ^,^(xi)^(x 2 )fi >, 

i 

A { +\x 2 -x 3 ) = < n,(P j {x 2 )xj{x z )n >, 
j 

A { 2 + \x 3 -x!) = n <n ' Mxs)Xk(xi)n > (5.1.10) 



the distributions appearing in the analysis of the second order perturbation theory i.e the lists ( 4.1.9 ) 



(4.1.12). They appear with these three combinations of arguments from various Wick contractions 



in the preceding formulas for the commutators. Here the fields <ft(xi), ip(x 2 ), x( x 3)) are factors 
in the Wick monomials of T{x\), T{x 2 ), T{x^) respectively. If Fermi fields are present one has to 
take into account the signs induced by the permutation of the non-commuting factors in defining the 
associated distributions A(~). 

We have to investigate the types of numerical distributions with causal support which can appear 
from the computation of the four commutators. These distributions will depend only of two variables 
£i = xi — x 3 , £2 = x 2 — x 3 due to translation invariance. It convenient to use again a graph theory 
terminology. We define a super-line to be the assemble of lines of a Feynman graph connecting two 
vertices. Then the notions of super-tree and super-loop are obvious and we have only such types of 
graphs. We give the generic form of the distributions associated to them. 

(a) First we obtain some distributions containing a factor 5 from commutators containing a factor 
5(x — y)L(x) or 5(x — y)L fJ, (x). In this case we obtain distributions of the type 

d(A)(x 1 ,x 2 ;x 3 ) = S(xi - x 2 )A(x 2 - x 3 ) (5.1.11) 

and other permutations of the variables. Here the distribution A is one from the lists ( gXp - fljllg ). 

(b) Next, from super-tree graphs we obtain three types of distribution. 

(bl) There exists a super-line between x\ and x 3 and a super-line between x 2 and x 3 . In this case 
one obtains distributions of the form: 



d(3)(xi,x 2 ;x 3 ) = A^ +) (x 2 -x 3 )A|, \x 3 -x 1 )-A[ \x 2 - x 3 )A 2 +) (x 3 - x\) 
+A%(x 3 - x 1 )Ai(x 2 - x 3 ) - Af(x 2 - x 3 )A 2 (x 3 - xi) 



(5.1.12) 



The causal support of this type of distribution can be checked if one derives a alternative for- 
mulae. If: 

^ (ret) (x 1 ,x 2 ;x 3 ) = Af ( ^ ) (x 3 -x 1 )Af t,(ret) (x 2 -x 3 ) (5.1.13) 



then we have from ( f4.1.46| ): 

d( 3 )(xi,x 2 ;x 3 ) = d^(x 1 ,x 2 ;x 3 ) - d r $(x 1 ,x 2 ;x 3 ). 



(5.1.14) 



Moreover, if one uses the expression of the third-order chronological product (3.1.22) one can prove 
that the distribution of this type is producing the Feynman propagator 



df 3) {xi,x 2 ;x 3 ) = A%(x 3 - xi)Af (x 2 - x 3 ) 



(5.1.15) 
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(b2) There exists a super-line between x\ and x 2 and a super- line between x\ and x 3 . In this case 
one obtains distributions of the form: 



d(i)(x 1 ,x 2 ;x 3 ) 



A 



(+), 



x s 
, AF 



xi)A { 3 \ Xl 



x 2 )-A { 2 \x 3 -x 1 )A^'(x 1 -x 2 ) 
+A3 * {x\ - x 2 )A 2 (x 3 - xi) - Af (x 3 - xi)A 3 (xi - x 2 ) 



Mi 



(5.1.16) 



The causal support of this type of distribution can be also checked if one derives the alternative 
formulas. We define: 



,adv(ret) / \ \ ret(adv) , « . ret(adv) , \ 

\x\,x 2 \x-i) = A 3 / (xi-x 2 )A 1 (x 3 -xi) 



(5.1.17) 



and we have as before: 

(aJi»a?2; x 3 ) = dffi{xi,x 2 ;x 3 ) - d r ^(xi,x 2 ;x 3 ). (5.1.18) 

If one uses the expression of the third-order chronological product ( 3,1.22| ) one can prove that the 
distribution of this type is producing the Feynman propagator 

df 1 - ) (xx,x 2 ;x 3 ) = Af (xi - x 2 )Af(x 3 - xi) 



(5.1.19) 



(b3) There exists a super-line between x\ and x 2 and a super- line between x 2 and x 3 . In this case 
one obtains distributions d( 2 )(xi, x 2 ; x 3 ) of the same form as in case (b2) if one makes xi <-> x 2 . 

We will denote the distributions associated to super-tree graphs by d^)(A , A')(xi, x 2 ; x 3 ) indi- 
cating explicitly the distributions in one variable A, A' from the lists ( 4.1.9f )-( [4.1.12| ) involved in 
the construction. One can verify that if the orders of singularity of these distributions are u and a/ 
respectively, then: 

w(d (i )(A, A')) = 4 + u + u/. (5.1.20) 

(c) We consider now graphs with a purely Bosonic super-loop. One obtains the following type of 
distribution; 



d(i 23 )(xi,x 2 ;x 3 ) = Af F (xi - x 2 ) A ( \ r \x 2 - x 3 )A { 2 \x 3 -xi)-A\ \x 2 - x 3 )A < 2 + \x 3 - xi) 

+Af (x 3 ) - xi ) [a 3 +) (xi - x 2 )aS _) (x 2 - x 3 ) - A^ } (xi - x 2 )AS +) (x 2 - x 3 ) 

+Af (x 2 - x 3 ) [a[ + \x 3 - xi)A^(xi - x 2 ) - A^ } (x 3 - xi)A< +) (zi - x 2 )] ; (5.1.21) 

for a Fermionic super-loop a overall —1 sign appears. 

The causal support property can be checked by deriving two alternative formulae : 



d(i 23 )(xi,x 2 ;x 3 ) 



x 2 ) A^ +) (x 2 - x 3 )A 2 ; (x 3 - xi 



(-)# 



A\ \x 2 



x 3 )A { 2 + \x 3 



+Af"(x 3 ) - Xl ) \Af\ Xl - x 2 )A ( -\x 2 - x 3 ) - A < f\x 1 - x 2 )A[ +) (x 2 - x 3 ) 



+Af v (x 2 -x 3 ) A^(x 3 



.(+), 



xi)A 3 } (xi 



-A^(xi 



x 2 ) aJ +) (x 2 -x 3 )A 



(-) 
2 



- x 2 ) 
(x 3 - 



+A r 2 et (x 3 ) - Xx) 

,(+), 



4 +) (x 1 -x 2 )a;- ) (x 2 

.(-) 



-A 

Xi) 

x 3 ) 



(-), 



^3 



X1)A( +) (X! 



x 2 . 



A^ } (x 2 
A^\ X1 



+A r f{x 2 -x 3 ) AY > {x 3 -x 1 )A\ s {x 1 -x 2 )-A ( -\x 3 -x l )Af>{x l -x 2 ) . (5.1.22) 



x 3 )A^ +) (x 3 

x 2 )a[ +) (x 2 
.Mi 



xx, 
X3) 
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We denote suggestively this type of distributions by dn 23 ) (Aj., A2, A 3 )(a;i, x 2 \ x 3 ) where Aj, i = 
1, 2, 3 are distributions from the lists ( 4.1.9| )- (|4.1,12| ) and we have concerning the order of singularity: 

w(d (123) (A 1) A 2 ,A 3 )) = 4 + ^a;(A i ). (5.1.23) 

i 

We say now something about the generic momentum space structure of such a distribution. First 
one has to obtain from the explicit formulas for the distributions A in one variable that in all cases: 

A i {± \p)^e(± Po )f i (p 2 ) (5.1.24) 

with supp(fi) C {p 2 > A 2 } for some parameters with mass significance Xi > 0, i = 1,2,3. We 
consider now the Taylor transform of A( 123 )(£i, £2) and we use the notation K = k\ + k 2 ; the generic 
structure is: 

A ( i23)(&i, h) = 0(kj - (A 2 + A 3 ) 2 )<?i + 0(k 2 2 - (A 3 + Xi) 2 )g 2 + 9(K 2 - (X 1 + A 2 ) 2 ) 53 - (5.1.25) 

It follows that if at least two of the masses Xi > 0, i = 1,2,3 are strictly positive, then (0,0) 
supp(A( 12 3)(/ci, k 2 ). This observation is useful because for causal distributions with such support 
property in momentum space one can use the so-called central formula for causal decomposition of 
distributions |23| ] . If the conditions of validity of the central formula are not meet we will have to use 
a regularization procedure. 

(ii) We investigate the possible Ward identities and obstructions to causal splitting. First we 
consider the case (b). We illustrate this case on the the distribution 

d" = d (3) (0"£> m ,A) (5.1.26) 

where A is arbitrary. The other cases can be treated similarly. First we derive the Ward identity: 

= -5( Xl - x 3 )A(x 2 - x 3 ) + m 2 D^( Xl - x 3 )A(x 2 - x 3 ). (5.1.27) 

Using the formula for the causal splitting ( |5.1,13| ) one can see that the preceding identity is 
preserved by the operation of distribution splitting. 

Next, we consider case (c). We have to study separately the case when the superdoop contains at 
most one Dirac line and the case when we have three Dirac lines. We illustrate the first case one the 
distribution 

d» = d {123) {d»D mi ,D m2 ,D m3 ), mi>0; (5.1.28) 
the other cases can be treated similarly. The Ward identity is in this case: 

= -8(x 2 - x 3 )D m2)mi {x 3 - xi) + • • • (5.1.29) 

where by ... we mean contributions with the order of singularity strictly smaller than 0. One computes 
immediately that both hand sides have the order of singularity equal to 1. If we have m 2 + m 3 > 
then we can apply the central decomposition formula and obtain no anomaly. In the opposite case, 
we use standard regularization procedure (|5.1.lj ) of the distributions appearing in the lists ( 4.1.9 )- 
(4.1.12) presented at the beginning of this Subsection. The decomposition (5.1.1) induces a similar 
decomposition for the distributions of the type d^y. 

d {{] = d° + d reg (5.1.30) 



28 



where 

uj(d°) = u(d), Lo(d re 9) = uj(d) - 2 (5.1.31) 

and the support properties of d° in the momentum space are more convenient: (0,0) supp(d^). 

If we apply this decomposition to the distributions cP and d we get two Ward identities, one for 
each piece. The first one can be split causally without anomalies using the central decomposition 
formula. For the second identity we note that both hand sides have order of singularity strictly lower 
than — 1 so this relation can be also split causally without anomalies as explained at the end of Section 

In this way we can obtain a anomaly-free decomposition of the Ward identity we have started with. 
One has to check case by case this argument for all the other types of distributions of type (c) without 
Dirac loops. 

A very important observation is that the preceding argument is not valid for distributions associ- 
ated to super-loops containing three Dirac lines. The reason is that one is lead to the computation of 



some traces. To be more specific the relevant terms from the first commutator ( 5.1.9 ) are: 



Df(x 1 ,x 2 ;x 3 ) = d^(x 1 ,x 2 ;x 3 ) : u a (xi)A bu (x 2 )A cp {x 3 ) : 
+ d abc( x ±> x 2' x 3) : u a {xi)^ b (x 2 )A cv (x 3 ) : + (x 2 <-> x 3 ) 
+d^ bc (x 1 ,x 2 ;x 3 ) : u a (xi)$ b (x 2 )$ c (x 3 ) : 
+d%(xi,x 2 ;x 3 )u a (x 1 ) H 



(5.1.32) 



where by • • • we mean the terms which cannot produce anomalies. Let us note that all these terms 
are obtain from Wick contractions of the pieces of the interaction Lagrangian of the type ( p. 2. 7 ). 

The distributions appearing in this formula are sums of distributions of the type d( 12 3) because of 
the traces. But in this case, the trace operation can annihilate the most singular term and instead of 



( 5.1.23 ) we might have: 



uj{d) < 4 + ^u;(Ai). 



(5.1.33) 



It follows that these distributions must be studied separately and some explicit computation are 
required. 

(iii) All the distributions appearing in the formula (|5 . 1 . 32 ) have 8 contributions corresponding to 
the decomposition of the three currents involved in ( |2.2.7| ) into the vectorial and axial components. If 
we compute the contribution corresponding to three vectorial factors, then the others can be obtained 
by simple substitutions. Let us consider this pure vector contribution to d^(xi,x 2 ;x 3 ); the explicit 
expression is: 



a abc:VVV 



{t a ) A B(tb)Bc{tc)cAd^% } AMB + {t a ) AB {t c ) B c{tb)cAd V M P c \ V M B ,M A (5.1.34) 
where we have defined for arbitrary masses M\,M 2 ,M 3 the following fundamental distribution: 

^Mi P ,M2,M 3 ( x l' x 2! ^3) = 

Tr{S$f 3 On - x 2 )Y \s { m ) 1 (x 2 - x 3 )YS[f 2 (x 3 - Xl ) - S$(X2 - x^S^ix, - Xl j\ ^ 

x 2 ) 7 



x sh p S { M ] 2 (x 3 - x 1 )^sif 3 (x 1 - x 2 ) - Sif 2 (x 3 - x l )YS { M l{x 1 



+ S F M2 {x 3 - Xl )^ S^(x 1 -x 2 )YS^>(x 2 -x 3 )-S^>(x 1 -x 2 )YS^(x 2 -x 3 ) Y} (5.1.35) 



:(+), 



?(-)/ 



which is similar to (5.1.21); compare also to formula (5.3.11) from |23| ]. It also has causal support: 
one can obtain quite easily alternative expressions having the structure ( |5.1,22| ). 
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The entire vectorial contribution is now obtained if we add the contributions following from 
^abc-VVV ^ we P er f° rm the following simple transforms: 

t -<, t b ^t' b , 7 M ^7 M 75, Y^Yls, (5.1.36) 

and the other two similar possibilities. Using the formula ( 4.1.13j ) we obtain the following form for 
the pure vectorial part: 

d abc;V = ( t a)AB(tb)Bc{tc)cAd^JJ AMB + (t a ) AB^c) BC^b)cA^M^M B M A 

+ {t' a )AB{t' b )Bc{tc)cAd^McM A ,MB + (4)AB(i C )Bc(4)cA^M c ^Mj3,Ai jl 
+ {ta)AB{A)Bc{t'c)cAd A ^ P J\^} MAMB + {t a )AB{t' c )Bc(t b )cAd^^ MB)MA 

HQAB{t h ) BC (t' c )cAdZfJ A ^ MB + (t' a )AB{t' c )Bc{tb)cAdZ P c {V J B ,- MA (5-1-37) 

One notices that the vectorial part of is express only in terms of the distribution of the type 

a M 1 ,M2,Ms' 

By similar transforms one can obtain the pure axial part. One has in defines in analogy to ( ^.1.34| ) 
the distribution d^M^ Mg by inserting a factor 75 : 

^mI,mM^ ^3) = Tr 75 {---} (5.1.38) 



where by {• • •} we mean the same paranthesis as in (5.1.35). The pure axial contribution to a\£ is 
similar to ( 5.1.37| ) . The only relevant thing is that it is expressed only in terms of the new distribution 
^mTm2 M3' ^ f° nows that the distribution can be expressed in terms of two independent 
distributions: d^^ M and d^Ma Mg . One can prove quite easily that the order of singularities are 

^ZSm\) = 1- (5-1-39) 
Let us note in passing that the asymptotic behaviour of the distribution 

dT = dT,v + dT,A (5-1.40) 

is given by: 

d^~V abc d^ )+ A abc d^ (5.1.41) 

where the axial tensor A abc is given by the expression ( |1.0.1| ) from the Introduction, the vector tensor 
is given by a similar expression: 

V abc = Tr(t+[t+, #] + t-[tfc,t-]) = fbcdTr{t+t+ + t-q) (5.1.42) 

and 

a (V)(A) = a 0,0,0 ■ {0.1.16) 

A similar investigation can be performed for the other distributions appearing in the formula 
(5.1.32). The distribution d^ bc can be expressed in terms of two independent distributions: d^^ 2 Ms 

and dj^^-j m 3 which can De obtained from 4/f m] M3 anc ^ ^MimI m 3 making 7 P — > 1. The order of 
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singularity of these distributions is lower than naive power counting indicates. They can be written 
as follows: 

1=1 

with 

U ^(V)(A) )=0 _ (5L45) 

Analogously, the distribution d^ bc can be expressed in terms of two independent distributions: 

d M^M 2 ,M 3 and d MuM 2 ,M s which can be obtained from d Mi,M 2 ,M 3 and d Mi,M 2 ,M 3 making Y -» 1. We 
also have: 

^4Vm% 3 ) = i- (5-1.46) 
Finally, the distribution da can be expressed in terms of two independent distributions: e^-Mi m 2 m 3 

and d m^M u M2,M3 which can be obtained from d^M^Ma and d MuM 2 ,M 3 b y making Smi £ m ,Mi- We 
have in this case: 

"(<S2r 2) M3) = 3- (5.1.47) 
We also need the distributions M ^ M ^ and df^) Mi M2 Ma which can be obtained from d^p M ^ M ^ 

and d^ jvli M 2 M 3 respectively by making 7^ — » 1. In this case we have a structure similar to ( |5.1.44j ) 
and a similar result for the order of singularity. 

We can easily see that all distributions d"'^ are completely antisymmetric in the Lorentz indices 
due to traces involving a 75. matrix. It is not difficult to prove that one can impose a supplementary 
condition on the causal splitting procedure, namely the preservation of this symmetry property. 

The distributions appearing into the third commutator from ( |5.1.9| ) can be obtained from the 
preceding ones by making the substitution x\ «-> x 2 and this doubles the value of the possible anomalies 
coming from the first commutator. 

The distributions appearing in to second commutator from fl5.1,9|) can be obtained from the pre- 
ceding ones by more subtle transforms. For the case (b) and the case (c) without Dirac loops we have 
the same list of distributions and there are no anomalies. For the case (c) with Dirac loops we have 
to consider: 

a M 1 ,M 2 ,M3 a M u M 2 ,M 3 
( ^M 1 ,M 2 ,M 3 ( Xl ' X 25 £3) — > /M 2) M 3 ,Mi( a ' 1 ' X 2' X3 ) = d M 2 ,MsMi ( X2 ' X3 ' Xl ) 
dlhM^Ms^ 1 ^ 2 '^^ fth,M 2 Mz^ x ^ x ^ ~ d M 3 ,Mx,M^ x ^ x ^ x ^- (5.1.48) 

(iv) Now we can give the list of Ward identities verified by these distributions. Using Dirac equation 
for the propagators we get: 



l ^ d M*Zl,MS Xl ' X2 '> X ^ = ( Ml " M 3)fMuM 2 ,M^ Xl ' X2 ^ X ^ 

iS(xi - x 2 )P^ 2M3 (x 1 - x 3 ) - i5{x l - x 3 )P^ uM2 (x 2 - x 3 ) (5.1.49) 



+iS(xt - x 3 )P^ iM2 {x 2 - x 3 ) + i5{x 2 - x 3 )P^ 3;Ml (xi - x 3 ) (5.1.50) 
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^ d Mi,M 2 ,M 3 ( X l' X 2;X3) = (Mi ~ M 3 )d^ M3 Ml (x 2 ,X 3 -Xl) 

-iS(xi - x 2 )Pm 3: m 2 {x 1 - x 3 ) - i5{xi - x 3 )Pl h M2 (x 2 - x 3 ) (5.1.51) 

^/S.Ms^.^i^) = (M 3 - M 2 )4^ 3iMi (x2,x 3 ;xi) 

+iS(xi - x 3 )Pl h M2 (x 2 - x 3 ) + i5(x 2 - x 3 )Pm iM3 {xi - x 3 ) (5.1.52) 

i ^ d m^,M 3 i x ^ x 2\x 3 ) = (Mi - M 3 )d ( ^l M3 Mi (x 1 ,x 2 ;x 3 ) 

+i5(x 1 - x 3 )P Ml ,M 2 {x 2 - x 3 ) + iS(x 2 - x 3 )P Mi ,m 3 (xi - x 3 ) (5.1.53) 



d 

i 



g^^S.Ma.Ma^l'^;^) = (Mi - M 3 )d^M 2 ,M 3 ,MS Xl > X2 ' X ») 

+i5{xi - x 3 )P m - Ml ,M 2 (x 2 - x 3 ) +i5{x 2 - x 3 )P m - Ml ,M z {xi - x 3 ). (5.1.54) 

The Ward identities for the axial distributions present a notable difference. Because of the trace 
operation, the delta terms disappear. Using also formula ( f4.1.13| ) we get: 

i ^ d Mi!&,M 3 ^ Xl ^ X2]X ^ = ( Ml + M 3)/m 1 ( i m 2i m 3 ( x 1' x ^ &s) (5.1.55) 

i ^ d MiM 2 ,M3^ Xl ' X2 ' X ^ = ( Ms + M 2)4/J A M 2 ,M^ XU X2! X ^ (5.1.56) 

i 7^ d Mi A M 2 ,M3 < ^ Xl ' X ^ X ^ = ( M i + m 3)4/^ ) M3,m 1 ( x 2,x 3 ;xi) (5.1.57) 

i^fhhShM^ 1 ^ 2 '^^ = ( M3 + M *)tf!&tz&S x 2> x & x l) (5.1.58) 



3 



I- 



~ d MuM2,MS Xl ' X2 '^ X ^ = (5.1.59) 



I- 



dx± 

9 



> dx » d m;M u M 2 ,M3( X ^ X2 '> X z) = °" (5.1.60) 

The causal splitting of these two type of Ward identities is sensibly different. Let us first consider 
only the first six equations (the vectorial Ward identities). Because of the delta terms in the right 
hand sides, we have the same order of singularity for both sides in all vectorial Ward identities , if the 
conditions of application of the central splitting formula are meet we get no anomalies. If some of the 
masses are null, one has to use a regularization procedure like for the case (b). More precisely one 
can prove that the decomposition ( 5.1.1] ) induces a similar decomposition for the distributions of the 



type d(i23) : 

d (123) = ^(123) + d (123) (5.1.61) 

where 

W Kl23)) = W ( d (123)), ^( d (lf 3 )) = W ( d (123)) - 4 (5.1.62) 



32 



and the support properties of d9 123 ^ in the momentum space are more convenient: (0, 0) supp(d) 12 3)) • 
We turn now to the last six equations (the axial Ward identities). The previous argument is still 
valid for the last two of them. We consider some generic anomalies P-f, Pg obtained after the causal 
splitting of the identities ( |5.1.57| ) and ( 5.1.58 ). If we differentiate the corresponding equations with 
respect to x v ^ and x\ respectively, we obtain the consistency equations 

Vi^O, A P 3=0 (5.1.63) 



dx\ 1 ' dx\ 

and this leads to P? = 0, i = 1, 3. 

The Ward identities ( 5.1.55 ) and ( 5.1.56| ) can produce anomalies of the type 



2 



P"P(X) = const E vpa ? - 5{X) (5.1.64) 

dxfdxi^ 



for some positive constant const. The explicit expression of this constant can be computed as in [23 
Section 5.3. One cannot eliminate such type of anomalies from both equations by redefinitions. The 
resulting anomaly is then: 

A(X) = const. A ahc e v ^ - 5(X) : u a {x 1 )A bv {x 2 )A cp (x 3 ) : (5.1.65) 

dxf 8x2 

where: 

A abc = 2 Tr(t' a {t' b ,t' c } + t a {t b ,t' c } + t' a {t b ,t c } + t a {t' b ,t' c }). (5.1.66) 

Performing some redefinitions of the expressions Af(X) ("integration by parts") we can reexpress 
this axial anomaly in the following form: 

Aabbj(X) = const. A abc e pvpa 5{X) : u a {x{)F^ { Xl )F^{x z ) : (5.1.67) 

and we can also show that the tensor depending only on the group indices A abc is in fact given by the 
formula ( |1.0.1| ) from the introduction. The anomaly Aabbj is a cocycle 

d Q A ABBJ = (5.1.68) 

but it is not a coboundary; so disappears iff we have the condition A abc = i.e. the well-known 
condition ( |5,1.6j ) from the statement. 

(v) We still have to investigate the possible anomalies coming from the delta terms i.e from distri- 
butions associated to graphs of type (a). We present here briefly the analysis of these terms. One can 
compute the commutators and select the terms which will lead, in principle, to an anomaly. We get: 

d 

[Tf (x),L{y)] = fabcfdcffdgh-Q^rDm c (x - y) ■ u a (x)A bu (x)A fX (y)Ag(y)A^(y) : 

d 

y) : u a {x)A bp {x)A p h {y)<5> e {y)<5> g {y) : 
+V'a b Jd b ff' dgh ^-D mt {x-y) : <$> a (x)u c {x)A fp {y)AP h {y)<5> g {y) : 
+if ahc g' b f 9 h-^D ml {x-y) : * a (x)u e (x)* f (y)%(y)* h (y) : + ■•• (5.1.69) 
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By • • • we mean the rest of the commutator which cannot produce anomalies Now, as in [16] and 
[17 1 we get from this commutator a possible anomaly: 



where 



A( Xl ,x 2 , x 3 ) = 5(X) L4i(xi) + A 2 {x x ) + A 3 (xi) + A^)] 

Ai(x) = ifabcfdcffdgh ■ u a (x)A bu (x)Afx{x)A u g (x)A^(x) : 
A 2 {x) = -2if abc f' dec f' dgh : u a (x)A bp (x)A p h (x)$ e (x)® g (x) : 
A 3 {x) = 2if abc f' dbf f dgh : <Z> a {x)u c (x)A fp {x)A p h (x)<I> g (x) : 

A±(x) = i [f' abc g' b f g h + f'fbcd'bagh + f'gbcd'bafh + f'hbcd'bafg] 

: <f> a (x)u c (x)$f(y)<f>g(x)<f> h (x) : 



(5.1.70) 

(5.1.71) 
(5.1.72) 
(5.1.73) 

(5.1.74) 



The results are: 



In [11 1 it is proved that the A\ = due to the Jacobi identity. 



• One can also show, using the identity (2.1.19) that A 2 + A3 = 0. 

• If we try to write the anomaly A4 as a coboundary cIqL(x) we should take 

L{x) = g' acfgh : ®a(x)$c(x)®f(x)$ g (x)$ h (x) : 
which is forbidden by the assumption that ( p.l.lCj ) is fulfilled. 
So we get the second restriction from the statement. I 



(5.1.75) 



Remark 5.2 Recently /^/ a new method was proposed to solve problems of consistency as the ones 
appearing in our paper. Instead of imposing a factorisation condition of the type ( 2.1.1(\ ) (or its 



"infinitesimal" version (2.1.11)) one imposes a quantum analogue of the Noether conservation law of 
a certain current. Presumably, this starting points are equivalent and they should lead to the same sets 
of consistency conditions. This point deserves further investigations. However, one should compare 



carefully the relation ( {.2.t ) expressing the conservation law of the BRST current ( and equivalent 
to the formal adiabatic limit condition) to the relation (4-5) of [TQ] expression the quantum Noether 
postulate. 

5.2 The Standard Model 



We remind the notations from [17]. The Lie algebra is in this case su{2) x u{l) and the standard basis 
X a , a = 0, 1, 2, 3 has the usual commutation relations 

[X a ,X b ] =e abc X c , a, b = 1,2, 3, [X ,X a }=0, a =1,2,3. (5.2.1) 

In the new basis Y a , a = 0, 1, 2, 3 defined by 

Y a = gX a , a = 1,2, Y 3 = -g cosO X 3 + g' sin 9X , Y = -g sinO X 3 - g' cos 6X . (5.2.2) 
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(here the angle 9, determined by the condition cos 9 > is the Weinberg angle and the constants g 
and g' are real with g > 0) the structure constants are: 

/210 = 9 sin 9, /321 = g cos 9, f 310 = 0, / 32 o = (5.2.3) 

and the rest of the constants are determined by antisymmetry. The choice of the masses is: 

m = 0, m a ^0, a =1,2, 3 (5.2.4) 

(the particles created by Aq being the photons and the particles created by Aa , a = 1,2,3 the heavy 
Bosons). 

In f\M we have found out the following result: 



Theorem 5.3 In the standard model, the following relations are true: 
(a) the masses of the heavy Bosons are constrained by: 



mi = ni2 = m^cos 9; (5.2.5) 



(b ) the constants f' abc are completely determined by the antisymmetry property ( 2.1.11 ) and: 

26 
1 

(5.2.6) 



f _ f > _ 6 9 ff _ e g , , _ , _g , _ cos 29 

Jon — J022 — J033 — 2~cos~9~' _ ^ Sm _ _ _ 2' 3 _ _ ^ 2cos 9' 



the rest of them being zero. Here e can take the values + or —. 
(c) the constants f " abc are (partially) determined by: 

fabc = (a,b,c = 1, 2, 3), /001 = /002 = /003 = /012 = /023 = /031 = 0' 

/on = /022 = /033 = i2mi ( m Q^ 2 ' (5.2.7) 

Moreover, one can fix e = +. 

Remark 5.4 In J]D/ a dual point of view is followed: one gives the masses of the heavy Bosons 
mi = 1712 7^ m 3 an d determines that the gauge algebra must be su(2) x it(l). 

We consider the minimal standard model containing only one generation of Dirac particles. In this 
case one takes in the generic formalism from the preceding Section: N = 2 and 

:)• <«■»> 

The components ^2 (resp. tpx) correspond to the electron (resp. the electronic neutrino) and m e 
is the electron mass. Remark that the neutrino mass is considered null. 
The choice for the representations is the following one: 

11 
H = -jp^x, 4 = 2 (~9 C0S da z + d' sin 91) , 

t+ = X -go2, t+ = - - (gsin 9a 3 + g 1 cos 91) (5.2.9) 

and 

t+ = t+ = 0, 4 = y sin 9, t£ = -y cos 9; (5.2.10) 

here are the Pauli matrices. The representation property (4.1.5) is fulfilled for any matrix y. 
However, we have the following elementary result: 
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Proposition 5.5 The interaction between the Dime field of the electron tp2 and the electromagnetic 
field Aq has the usual form 

e : ^27/^2 : 

(here e is the electron charge) iff 

9=-?~Z, 9' = (5-2.11) 

sin v cos i) 

and 

y = ~«/(l - a 3 ) = ^—g'M. (5.2.12) 
2 2m e 

Next, we have 

Proposition 5.6 The expressions for the matrices s^ are: 

m e / 0\ , im e / 0\ . m e / 0\ . im e / 
1 1 sT = I „ , sj = . „ , s ' 



~2miV0 1J' 1 ~2m 1 \l Oj' 2 ~ 2m x \ 1 OJ' 3 2mi \0 1 

(5.2.13) 

Proof: 



One uses the relations ( 2.2. 18| ) for a = 1,2,3 and obtains the expressions for sj, a = 1,2,3. 
Next, we use the relation ( 4.1.6|) , more precisely: 

ta4 ~ 4*i = ifL*t> a = 1,2, 3. (5.2.14) 
This equation gives immediately the expression for Sq. M 

The expression of the Higgs potential is obtained as in [||. One can check that in this way 
the usual standard model is obtained. 



5.3 Regularization and Anomalies 

We have succeeded to give a complete analysis of the possible anomalies appearing in the standard 
model up to the order three of the perturbation theory. One would want to generalise this analysis to 
all orders of the perturbation theory. It is possible that one can use the same type of combinatorial 
argument, namely one considers possible distributions appearing in the commutators D(X) of order n 
and observes that only the super-loop graphs with Dirac lines can produce anomalies. Then it is quite 
possible that in higher orders the orders of singularity are sufficiently lower to make possible a causal 
splitting of the Ward identities without anomalies. This seems to be indicated by the traditional 
argument from the literature Q , , [|^] . 
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